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Abstract 

m 

CN ■ We present recent results on scalar light mesons based on Dalitz plot anal- 

yses of charm decays from Fermilab experiment E791. Scalar mesons are 
| found to have large contributions to the decays studied, D + — > K~it + it + and 

D + , Df — > 7r _ 7r + 7r + . From the first decay, we find good evidence for the exis- 
tence of the light and broad k meson and we measure its mass and width. We 



O 



OO ! find strong evidence for the cr(500) meson from D + — > 7r~7r + 7r + decay and mea- 

sure its mass and width. We also present the results obtained for the /o(980) 
parameters through the Df — > 7r~7r + 7r + decay. These results demonstrate the 



o 



importance of charm decays as a new environment for the study of light meson 



fl^ ; physics. 

& ■ 1 Introduction 

The decays of charm mesons are currently a new source of information for the study of 
light meson spectroscopy, with the advantages of having well defined initial state (the 
D meson, a 0~ state with defined mass). This new information is complementary to 
that from scattering experiments and can be particularly relevant to the understand- 
ing of the scalar sector, due to the well-known difficulties of scattering experiments 
to observe light and broad scalar resonances. 

Here we present an overview of the results we obtained analysing the decays 
D+ -> 7r-7r+7r+ 0, D+ -> it-it + it + § and D + -> A^7r+7r+ §, using data collected 
in 1991/92 by Fermilab experiment E791 from 500 GeV/c 7r - -nucleon interactions. 
For details see Q. 

For the D + 7r _ 7r + 7r + decays, we find that a model with only known 7T7T reso- 
nances plus a non-resonant (NR) channel is not able to describe the data adequately. 

1 To appear in proceedings of VIII International Workshop on Hadron Physics 2002 
Rio Grande do Sul, Brazil. World Scientific edition. 
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Thus we include a new amplitude in our fit function and find strong evidence for the 
presence of a light and broad scalar resonance, the <r(500). The channel involving 
this scalar meson is responsible for half of the decay rate. We measure the mass and 
the width of this scalar meson to be 478±i| ± 17 MeV/c 2 and 324j^ ± 21 MeV/c 2 , 
respectively. Other experiments have presented controversial evidence for this low- 
mass 7T7T resonance in partial wave analyses |5|, ||, [7[], with ambiguous interpretations 
of the characteristics of such particles[|], |§. 

We found a similar situation for the D + — > K~tt + tt + analysis. When we include 
all known Kit resonant channels plus a NR contribution, we find that this model 
is not able to describe the data. By including an extra scalar resonant state, with 
unconstrained mass and width, we obtain a fit which is substantially superior to that 
without this state. The values for its mass and width are found to be 797 ± 19 ± 42 
MeV/c 2 and 410 ± 43 ± 85 MeV/c 2 respectively. We refer to this state as the k. We 
also obtain new measurements for the mass and the width of the _K"q(1430) resonance. 

From our analysis of — > 7r~7r + 7r + decays we found that the isoscalar inter- 
mediate states are dominant in this decay as well. The largest contribution to this 
final state comes from the decay involving the scalar meson f Q (980) whose nature 
is a long-standing puzzle. It has been described as a qq state , a KK molecule, a 
glueball, and a multiquark state || |^]. We obtain new measurements for the /o(980) 
and / (1370) masses and widths. 

To obtain these results, we had to introduce a new approach in Dalitz plot analysis 
in order to extract the mass and width of the scalar resonances by using them as 
floating parameters in the fit. We begin this paper presenting the general method, 
applied in the D + — > K~tt + tt + Dalitz-plot analysis, then we discuss the — *> 
7r~7r + 7r + and D + — > tt~7t + tt + studies using the same procedure. 

2 The D + K-n + 7T + Dalitz-plot Analysis 

From the original 2 x 10 1 events collected by E791, and after reconstruction and 
selection criteria, we obtained the D + — > K~tt + tt + sample shown in Figure 0(a). The 
filled area represents the level of background; besides the combinatorial, the other 
main source of background comes from the reflection of the decay —> K~K + tt + 
(through K K + and </>7r + ). The cross-hatched region contains the events selected 
for the Dalitz-plot analysis. There are 15090 events in this sample, of which 6% are 
background. 

Figure 0(b) shows the Dalitz-plot for these events. The two axes are the squared 
invariant-mass combinations for K7r, and the plot is symmetrized with respect to the 
two identical pions. The plot presents a rich structure, where we can observe the clear 
bands from K (890)7r + , and an accumulation of events at the upper edge of the diag- 
onal, due to heavier resonances. To study the resonant substructure, we perform an 
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Figure 1: (a) The K 7r + 7r + invariant mass spectrum. The filled area is background; 
(b) Dalitz plot corresponding to the events in the cross-hatched area of (a). 

unbinned maximun-likelihood fit to the data, with probability distribution functions 
(PDF's) for both signal and background sources. In particular, for each candidate 
event, the signal PDF is written as the square of the total physical amplitude A and 
it is weighted for the acceptance across the Dalitz plot (obtained by Monte Carlo 
(MC)) and by the level of signal to background for each event, as given by the line 
shape of Figure fj](a). The background PDF's (levels and shapes) are fixed for the 
Dalitz-plot fit, according to MC and data studies. 

We begin describing our first approach to fit the data, which represents the conven- 
tional Dalitz-plot analysis including the known Kn resonant amplitudes (A n , n > 1), 
plus a constant non-resonant contribution. The signal amplitude, a Breit-Wigner 
parametrization times a Blatt-Weisskopf damping factor, is described in reference [[| . 

Using this model with well-known resonances (Model A), we find contributions 
from the following channels: the non-resonant, responsible for more than 90% of 
the decay rate, followed by Z*(1430)7r+, Z*(892)tt+, F*(1680)vr+ and F*(1430)vr + . 
The decay fractions and relative phases are shown in Table |I[ These values are in 
accordance with previous results from E691 |H| and E687 [pjj . We thus confirm a 
high non-resonant contribution according to this model, which is totally unusual in 
D decays. Besides, there is an important destructive interference pattern, since all 
fractions add up to 140 %. 

To evaluate the fit quality, we compute a x 2 from binned, two-dimensional distri- 
butions of data and decay model events. The \ 2 comes from the differences in the 
binned numbers of events between the data and the model (from a fast MC simu- 
lation). We obtain x 2 jv = 2.7 {y being the number of degrees of freedom), with a 
corresponding confidence level (CL) of 10 -11 . In Figure 0(a) we show the Kn low 
and high squared-mass projections for data (error bars) and model (solid line). The 
discrepancies are evident in the very low-mass region for m?(Kit)i 0W and near 2.5 



Decay 
Mode 


Model A: No re 
Fraction (%) Phase 


Model B: With re 
Fraction (%) Phase 


NR 

iT(892)7r+ 
^o(1430)tt+ 
K* 2 {U30)tt + 
K r *(1680)vr+ 


90.9 ± 2.6 0° (fixed) 

13.8 ±0.5 (54 ±2)° 
30.6 ±1.6 (54 ±2)° 
0.4 ±0.1 (33 ±8)° 
3.2 ±0.3 (66 ±3)° 


13.0 ±5.8 ±2.6 (349 ±14 ±8)° 
47.8 ± 12.1 ±3.7 (187±8±17)° 
12.3 ±1.0 ±0.9 0° (fixed) 
12.5 ±1.4 ±0.4 (48±7±10)° 
0.5 ±0.1 ±0.2 (306 ±8 ±6)° 
2.5 ±0.7 ±0.2 (28 ±13 ±15)° 



Table 1: Results without re (Model A) and with re (Model B). 



(GeV/c 2 ) 2 for m 2 (K7c)hi g h- These regions of disagreement are the same observed pre- 
viously by E687 |llj. We thus conclude that a model with the known Kn resonances, 
plus a non-resonant amplitude, is not able to describe the D + — > K~h + 'k + Dalitz 
plot satisfactorily. 

A similar pattern - bad fit quality with large NR fraction - is found in the analysis 
of the decay D + — > 7r~7r + 7r + when allowing only the established nir resonances [0 
. There we find that the inclusion of an extra scalar resonance improves the fit 
substantially, giving strong evidence for the a(500). See the section on D + — > 7t~7t + tt + 
below. Thus, we are led to try an extra scalar resonance in our fit model here. 

This second fit model, Model B, is constructed by the inclusion of an extra scalar 
state, with unconstrained mass and width. For consistency, the mass and width of 
the other scalar state, the _K"q(1430), are also free parameters of the fit. We adopt 
a better description for these scalar states by introducing gaussian-type form-factors 
[12] to take into account the finite size of the decaying mesons. Two extra floating 
parameters are the meson radii rjj and introduced for these meson sizes. 

Using this model, we obtain the values of 797 ± 19 ± 42 MeV/c 2 for the mass and 
410 ± 43 ± 85 MeV/c 2 for the width of the new scalar state (first error statistical, 
second error systematic), referred to here as the re. The values of mass and width 
obtained for the 7^(1430) are respectively 1459 ± 7 ± 6 MeV/c 2 and 175 ± 12 ± 12 
MeV/c 2 , appearing heavier and narrower than presented by the PDG pi. The decay 
fractions and relative phases for Model B, with statitical errors, are given in Table |l[ 
Compared to the results of Model A (without re), the non- resonant mode drops from 
over 90% to 13 ± 6%. The mr + state is now the dominant channel with decay fraction 
about 50%. The meson radii rjj and are found to be respectively 5.0 ± 0.5 GeV -1 
and 1.6 ± 1.3 GeV -1 , in agreement with values used by other groups [[1^, [14}[ . 

Moreover, the fit quality of Model B is substantially superior to that of Model A. 
The x 2 / v is now 0.73 with a CL of 95%. The very good agreement between the model 
and the data can be seen in the projections of Figure @(b). A number of studies were 
done to check these results. They are described in detail in reference B. 
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Figure 2: m 2 (K7Ti ow ) and m 2 (i^7rhigh) projections for data (error bars) and fast MC 
(solid line): (a) fit to Model A, without k, and (b) fit to Model B, with k. 



The Df 7Y-7T + 7T + Results 



s 



In Figure |3| we show the 7r~7r + 7r + invariant mass distribution after reconstruction 
and selection criteria [[l], [2[] for the sample collected by E791. Besides combinatorial 
background, reflections from the decays D + —> K~7t + tt + , D° — > K~n + (plus one 
extra 7r + track) and Df — > t]'ti + , v( — > p°(770)7 are all taken into account. The 
hatched regions in Figure |3| show the samples used for the Dalitz-plot analyses. There 
are 1686 and 937 candidate events for D + and Df respectively, with a signal to 
background ratio of about 2:1. The Dalitz plots for these events are shown in Figure |], 
the axes corresponding to the two 7r~7r + invariant-masses squared. 

For the Df — > 7r~7r + 7r + events in Figure |](a), the signal amplitude includes all 
channels with well-known dipion resonances ||: p°(770)7r + , / (980)7r + , / 2 (1270)7r + , 
/o(1370)7r + , p°(1450)7r + and the non-resonant, assumed constant across the Dalitz 
plot. 

For the /o(980) 7r+ amplitude, instead of a simple Breit-Wigner form, we use a 



coupled-channel Breit-Wigner function [15], 



BW fo{980) mL _ m 2 + zmo(rv + rV) 



.</ \ m'Ui - ml, T K = f- yfiqji - ml + + y/m^/A-m^) . (2) 
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M(n" 71* tt*) GeV/c 2 

Figure 3: The 7r~7r + 7r + invariant mass spectrum. The dashed line represent the total 
background. Events used for the Dalitz analyses are in the hatched areas. 

The — > 7r~7r + 7r + Dalitz plot is fit to obtain not only the decay fractions and 
phases of the possible sub-channels, but also the parameters of the /o(980) state, g n , 
gKi and mo, as well as the mass and width of the /o(1370). The other resonance 
masses and widths are taken from the PDG||. The resulting fractions and phases 
are given in Table [2|. The measured /b(980) parameters are mo = 977±3±2 MeV/c 2 , 
ftr = 0.09 ± 0.01 ± 0.01 and g K = 0.02 ± 0.04 ± 0.03. By fitting the Dalitz plot using 
a simple Breit-Wigner function, for the /o(980) we find m = 975 ± 3 MeV/c 2 and 
T = 44±2±2 MeV/c 2 , and the results for fractions and phases are indistinguishable. 

The confidence level of the fit for — > 7r~7r + 7r + is 35%. In Figure [5] we show 
the 7r~7r + mass-squared projections for data (points) and model (solid lines, from 
fast-MC). 



Decay Mode 


Fraction (%) Phase 


/o(980)tt+ 
NR 

P (770)tt+ 
/ 2 (1270)tt+ 
/o(1370)tt+ 
P (1450)tt + 


56.5 ±4.3 ±4.7 0° (fixed) 
0.5 ±1.4 ±1.7 (181 ±94 ±51)° 
5.8 ±2.3 ±3.7 (109 ±24 ±5)° 
19.7 ±3.3 ±0.6 (133 ±13 ±28)° 
32.4 ±7.7 ±1.9 (198 ±19 ±27)° 
4.4 ±2.1 ±0.2 (162 ±26 ±17)° 



Table 2: Dalitz fit results for Df — > n n + n 
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Figure 4: (a) The —> n Dalitz plot and (b) the D + — > 7T~ 

plot. Since there are two identical pions, the plots are symmetrized. 



7r + 7r + Dalitz 



4 The D 
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7T 7T + 7T + Results 



In a first approach, we try to fit the D + — ► 7r~7r + 7r + Dalitz plot of Figure f|(b) with 
the same amplitudes used for the Z?+ — > 7r~7r + 7r + analysis. Using this model, the 
non-resonant, the p°(1450)7r + , and the p°(770)7r + amplitudes are found to dominate, 
as shown in Table |3], and in agreement with previous reported analyses ]16, 17]] . 
However, this model does not describe the data satisfactorily, especially at low tt~tt + 
mass squared, as can be seen from Fig. ^(a). The x 2 / u obtained from the binned 
Dalitz plot for this model is 1.6, with a CL less than 10~ 5 . 

To investigate the possibility that another tt~tt + resonance contributes to the 
D + — ■> 7r~7r + 7r + decay, we add an extra scalar resonance amplitude to the signal 
PDF, with mass and width as floating parameters in the fit. 



Decay 
Mode 


Fit without ct(500)tt + 
Fraction (%) Phase 


Fit with a(500)7r + 
Fraction (%) Phase 


(t(500)7t + 
p°(770)7r+ 

NR 
/o(980)tt+ 
/ 2 (1270)tt+ 
/o(1370)tt+ 
P (1450)tt + 


20.8 ± 2.4 0° (fixed) 

38.6 ±9.7 (150 ±12)° 
7.4 ±1.4 (152 ±16)° 
6.3 ±1.9 (103 ±16)° 

10.7 ±3.1 (143 ±10)° 
22.6 ±3.7 (46 ±15)° 


46.3 ±9.0 ±2.1 (206 ±8 ±5)° 
33.6 ± 3.2 ± 2.2 0° (fixed) 
7.8 ±6.0 ±2.7 (57 ±20 ±6)° 

6.2 ±1.3 ±0.4 (165 ±11 ±3)° 

19.4 ±2.5 ±0.4 (57±8±3)° 

2.3 ±1.5 ±0.8 (105±18±1)° 
0.7 ±0.7 ±0.3 (319 ±39 ±11)° 



Table 3: Dalitz fit results for D + —> n 7r + 7r + 
systematics (only for fit with cr(500)7r + mode). 



First errors are statistical, second 



We find that this model improves our fit substantially. The mass and the width of 
the extra scalar state are found to be 478±|| ±17 MeV/c 2 and 3241^ ± 21 MeV/c 2 , 




Figure 5: s i2 and Si 3 (m% n ) projections for — > n data (dots) and our best 

fit (solid). The area under the dashed curve corresponds to background. 



respectively. Refering to this state as the a(500), we obtain that cr(500)7r + channel 
produces the largest decay fraction, as shown in Table [| The non-resonant amplitude, 
which is dominant in the model without cr(500)7r + , drops substantially. This model 
describes the data much better, as can be seen by the im mass squared projection in 
Fig. [|(b). The x 2 ' l v is now 0.9, with a corresponding confidence level of 91%. 




Figure 6: Su and S13 (rrv^) projections for D + — > -K~it + it + data (dots) and our 
best fit (solid) for models (a) without and (b) with cr(500)7r + amplitude. The dashed 
distributions corresponds to the expected background levels. 
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5 Conclusion 



From the data of the Fermilab E791 experiment, we studied the Dalitz plots of the 
decays D + — > K~7i + tt + , —> 7r~7r + 7r + and D + — > 7r~7r + 7r + . In these three final 
states, the scalar intermediate resonances were found to give the main contribution to 
the decay rates. We obtained strong evidence for the existence of <r(500) and k scalar 
mesons, measuring their masses and widths. We also obtained new measurements for 
masses and widths of the other scalars studied, /o(980), /o(1430) and i^o(1430). 

The results presented here show the potential of D meson decays for the study of 
light meson spectroscopy, in particular in the scalar sector. 
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